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^^ ■ Abstract 

,S^ • In this paper, we compute the natural density of the set of fc x rt integer matrices that 

j^ I can be extended to an invertible n x n matrix over the integers. As a corollary, we find 

the density of rectangular matrices with Hermite normal form [Okx(n-k) Ik] ■ Connections 
with Cesaro's Theorem on the density of coprime integers and Quillen-Suslin's Theorem 
are also presented. 

(N 

t^^ . Key Words: Natural density, unimodular matrices, Cesaro's Theorem, Quillen-Suslin's The- 

Qv I orem 

cn : Subject Classification: 15B36; 11C20 

o 

1 Introduction and Main Result 

Given a commutative ring R with 1, the notion of invertible nx n matrix is well defined, and 
^ ' can be characterized by the condition that the determinant of such a matrix is a unit in R. 

j^ ■ Given a kxn matrix A, the question of whether it can be completed by an {n — k)xn matrix 

into an n X n invertible matrix over R has raised several interesting problems in the past. For 
instance, the celebrated Quillen-Suslin Theorem [14^ I15j. previously known under the name 
Serre Conjecture, deals with the case when R = ¥[xi, ...,xz], with F a field and A; = 1. It can 
be shown that this case also contains in essence the general case 1 < k < n, see [18]. The 
theorem states that over this ring, the following three properties are equivalent |14 1 ll5 t [T8]: 

1. A can be completed into an n x n invertible matrix, 

2. there exists an n x A: matrix B such that AB = Ik, where Ik is the kxk identity matrix, 

3. the kxk minors of A have no common zeros. 

When the ring i? is a PID it is a direct consequence of the Smith normal form (see e.g. [TT]) 
that Conditions 1. and 2. are again equivalent and these conditions are equivalent to the fact 
that the gcd of the kxk minors is equal to 1. 
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In the sequel we will adopt the usual convention and call a k x n matrix A over some ring 
R unimodular as soon as A can be extended to n x n invertible matrix. 

Over the ring of integers various related results on unimodular matrices are known in the 
literature. E.g. Zhan |19j showed that any partial nxn matrix with n given entries not lying 
on the same row or column can be completed into a unimodular matrix. Another result is 
due to Fang {6j who showed that, if the diagonal of a square matrix is let free, then it can be 
completed into a unimodular one. 

Our focus in this paper will be on the "probability" that a random k x n integer matrix 
is unimodular. Related to our problem is a classical result due to Cesaro [H [21 [3] (see 
also |16l [T7] and the historical remarks below) which states that the "probability" that two 
randomly chosen integers are coprime is ^^^T ~ "^' where C denotes Riemann's zeta function. 
A re-statement of Cesaro's result is then: the probability that a random 1x2 integer matrix 
is unimodular is C(2)~^. 

In order to make the notion of probability precise we first remark that the uniform distri- 
bution over the set Z"^ has little meaning. For this reason researchers often use the concept 
of natural density when stating probability results in Z™ or more general infinite modules 
R"^. In the following we briefly explain this concept. Let S C Z™ be a set. Define the upper 
(respectively lower) natural density as 

l(5)=limsup '^'^J~i^'^["' , D(5)=liminf'^''[-^'^["' 



When both limits are equal one defines the natural density of the set S as: 

B(5) :=1(5) = D(S). (1) 

The following properties of natural density are readily verified: If S'^ denotes the comple- 
ment of S then D(5"^) = 1 — ©(S*), whenever S has a well-defined density. Similarly if {Si}i^i 
is a set of subsets of Z™ with well defined densities B(5i) and if S" = Ujg/S'j, then: 



KS)<Y.^iS,). (2) 



Readers interested in more background on the notion of natural densities of sets of integer 
matrices are referred to [8J, even though this paper attributes the result of Cesaro to Mertens, 
something which was probably triggered by a remark in |7j. 

In fact the exact fatherhood of the result of Cesaro appears to be inexactly described in 
several occasions in the literature. The story behind this historical misunderstanding seems 
to be the following. The problem of evaluating the probability that two random integers are 
coprime appears in a 1881 question raised by Cesaro [T|. Two years later, Sylvester [16] and 
Cesaro [2] independently publish their solutions. Interestingly, two proofs are presented in 
|16] : Sylvester's own argument is based on Farey series, and a more "probabilistic" argument 
of Franklin is also presented with his permission. In the footnote of a 1888 paper |17j . 
Sylvester publishes a similar proof, and mentions that Cesaro "claimed the prior publication " 
of the result. Sylvester's argument is based on Farey series, and the remark in |7] makes a 
connection between an earlier work of Mertens [12j of 1874 on the average value of Euler if 
function and Farey series, which probably triggered the remark in [8] but at no point in [7] 
the result is attributed to Mertens. If we want to associate the aforementioned probability 



with the average value of c/?, then it is legitimate to go back to an 1849 paper of Dirichlet 
[5] where the value -^ appears for the first time. The case of k coprime integers, /c > 2, 
is also presented for first time by Cesaro in 1884 [3]. The result is rediscovered in 1900 by 
D.N.Lehmer [9], apparently independently (see also [T3j). 

The major result of our paper is a matrix version of Cesaro's theorem: 

Proposition 1 Let 1 < k < n be integers. The natural density d^^n of k x n unimodular 
matrices with integer entries is given by 

dfc,n = ( n ^(■^■) 

vj=n—k+l 

In other words, the "probability" that a "random" k x n integer matrix can be extended 
into a matrix in GL„(Z) is given by (CC'^) • • • • • Cl'^ — k + 1))^ . Since Cl'^) converges rapidly 
towards 1, if fc = n — d, the above density converges rather fast to the limit dd with 



new 



0.43575707677... and ds = C(2) • di, dg = C(2) • C(3) • di 



When n is not too small, say n > 4, di is a good approximation of the proportion of integer 
matrices that can be completed by a row into an invertible matrix. We would like to point 
out that the proof of the above proposition is independent from the result of Cesaro and as 
such a new proof of his theorem is given when the above proposition is considered with k = 1. 
Proposition [1] above gives a probabilistic extension of the simplest case of the Quillen-Suslin 
Theorem. The concept of natural density does not extend naturally to the ring F[xi, . . . , x„] 
and thus the existence of a direct extension of our result to the general case is unclear. In 
the next section we will prove Proposition [TJ 

2 Proof of the Main Result 

Let us fix some notations for the rest of the article. Let l<t<A;<nbe integers. Given 
a k X n matrix A with integer entries, the determinants of the (") (j) t x t submatrices of A 
formed by the intersection of any subset of t column and t row vectors of A are called the 
i-minors of A. When t = k, they are called the full rank minors of A. The set of primes is 
denoted by P and when no confusion is possible, the set Z will be identified with Z . 

The strategy of the proof of the above proposition is to localize the computation of the 
density at every prime, and then lift the information up in order to extract the exact density 
over Z . The proof of the above proposition relies on the next lemmas. As noted before, the 
PID version of Quillen-Suslin's Theorem gives directly the following lemma: 

Lemma 2 A k x n matrix A with integer entries is unimodular if and only if the full-rank 
minors of A are coprime. 

Over a field F, the rank of a fc x n matrix, i.e., the maximal number of rows that are 
linearly independent over F, is equal to the largest integer t such that there exists a non-zero 



t X t minor. A classical computation shows that over the prime finite field F = Z/pZ, if Fp is 
the set of full rank k x n matrices, its cardinality satisfies \Fp\ = YljZoiP^'' ~ P')- See [10] for 
the details. This can be extended as follows: 

Lemma 3 Let S be a finite set of prime numbers. The density of kxn matrices with integer 
entries for which the gcd of the full rank minors are coprime to all primes in S is given by 

n n(i-^ 

j=n-k+lpeS ^ ^ 

Proof: Let us call Es C Z the set oi kxn matrices for which the gcd of the full rank minors 
are coprime to all primes in S. Let A^ = JliDGS^ ^^'^ recall the Chinese remainder theorem 
(Z/AZ) = ripes (^/p^)- Let B be an integer that will go to infinity in the sequel. Write 
B = qN + r, with g, r€N, 0<r<A and consider the map (j) obtained as the composition 
of maps 

[-gA,gA['="^ (Z/AZ)'=" -^ Yl (^/P^)''" 

p€S 

where the first map is the quotient modulo A and the second is the induced homomorphism 
given by the Chinese remainder theorem. Because of the above remark we have 

^{Esn[-qN,qNf^'') = YlFp. 

PG5 

The first map is a (2^) -to-1 map, i.e., each fiber contains exactly {2q) elements, and the 
second map is an isomorphism. Thus 

\Esf^[-qN,qN[^'^^\ = {2qf^■\{\Fp\. 

pes 

We have the disjoint union [-B,B[''"'= [-gA,gA['^'"U ([-S, B['^'"\[-gA,gA[^"). The differ- 
ence of hypercubes has a volume bounded by the volume of 2kn hyper-rectangles of side area 
^fcn-i ^j^j height r which gives < \[-B,B[^'^\[-qN,qN['^"'\ < 2knrB''^~^ and therefore we 
have 

\Es n [-B, Bf"'\ = \Es n [-(?Ar,gA['="| + p 

with 0< p< 2A;nr5^"-i. Thus 



l^sn [-5,5p| l^sn [-gA,gAp| {2qN) 



+ 



(25)fcn (2gA^)^" (25)^" {2B)f''' 

-w—r I F I / r \ kn 

pes ^ 



Finally, we have 



pes ^ pes ^ j=n-k+ipes ^ ^ -' 



D 



Extending the definition of Es, let us call Et C Z the set oi k x n matrices for which 
the gcd of the full rank minors are coprime with the first t primes 2,3, ... ,pt. Note that 
the sequence Et is a decreasing sequence of sets, i.e., Ei C Ej if i > j and ii E = PligN-E't, 
then Lemma [2] implies that E is the set of A: x n unimodular matrices with integer entries. 
In order to prove Proposition [H we have to prove that 0{E) exists and compute its value 
dfc,n = 9{E). Since we know B(£^t) for all t, it is tempting to prove the proposition by simply 
letting t going to infinity in the expression given by Lemma [3l but this is an invalid argument 
in general. Indeed the example of Et = [t,oo[ with E = % shows that it is possible to have 
a sequence of sets Ei C Ej '\i i > j with E = Ht^nEt and D{E) = 7^ 1 = linii^oo ^{Et) 
since IJ){Et) = 1, Vt. The next lemma describes how to avoid this pathological case. Let 
us recall that for any real sequences an,bn, liminfa„ + liminfi^ < liminf(a„ + 5„) and 
lim sup ( On + bn) < lim sup a„ + lim sup bn and lim sup — a„ = — lim inf a„ 



^n- 



Lemma 4 Let Et be a sequence of decreasing sets in Z"* such that D{Et) exists for all t and 
converges to d. Let E = Ht^nEt- If limt^oo^iEt \E) = 0, then D(E) exists and is equal to 
d. 

Proof: We will use the disjoint union Et = E U (Et \ E). Since 

\Etn[-B,Br\ = \E n[-B, Br\ + \{Et\ E) n[-B, Br\, 

we have 

\En[-B,B["'\ _ |^in[-S,S[™| -\{Et\E)n[-B,B["'\ 



(25)™ (2B)™ (25)™ 

Taking the lim inf and the lim sup, and since liniB^oo * AbV" ~ ^{Et), we have 

0{Et) - 0{Et \E)= B{E) < 0{E) = B(^j) - B{Et \ E). 
The result follows when t -^ cx) since < D(^i \E)< ^{Et \E) — ^ 0. D 

Proof of Proposition \^ We will use the previous lemma. The set Et\ E is the set of fc x n 
matrices A for which there exists a prime p with p > pt so that p divides the gcd of the full 
rank minors of A. For each prime p, let us define Hp to be the set oi k x n matrices whose 
gcd of the full rank minors is divisible by p. Then Et\E = Upyp^Hp. Note that Lemma [3] 
applied to 5 = {p} implies that Hp has a density equal to 



iHp)=i-o{H;)=i- n (i-^)- 

j=n~k+l ^ ^ ^ 



By induction on the number of factors, one readily verifies that for real numbers < x^ < 1, 
we have the inequality n?=o(-'- ~ ^i) > 1 ~ S?=o ^i' which applied to the above product gives 



7=n— fe+l 7=»i— fe+1 



j=n- 

Finally, we have 



1 2 

pj' p"-^(p — 1) p2 



\Et \ E) = B{Up^p,Hp) < Y, HHp) < E 4 

p>pt p>pt 



which shows that hnit^oo 0{Et \E) = since the last series is the tail of the convergent series 
X^pgp ~^- We can therefore apply the previous lemma and conclude that dfc^„ = D(i?) exists 
and is equal to 

j=n-k+lpeP ^ ^ j=n-k+l 

D 



3 Concluding Remarks and Further Results 

Proposition [T] does not cover the square case k = n. Indeed, the zeta function is not defined 
when n — /c + 1 = 1 since C, has a pole of order 1 at x = 1. The following result covers the 
square case. 

Lemma 5 The natural density of n x n unimodular matrices is d„,^„ = 0. 

Before we give a proof we remark that the product formula in Proposition [T] naturally has an 
extension to zero as li'nix-^i{C{^))~^ = 0- 

Proof: Each nxn matrix with n^ — 1 entries in the range [—B, B[ can be completed by at most 
two values in order for this matrix to be unimodular, due to the Lagrange expansion of the 
determinant, which must be ±1. As such there are at most 2(25)" ~^ unimodular matrices 
with entries in [—B, B[. The conclusion follows since dn,n < limB^oo 2(2i?)"("^^)/(2i?)" = 0. 

D 

The result of Proposition [T] can also be used in the determination of the natural density 
of k X n matrices whose Hermite normal form (HNF) is very simple. Recall that the HNF of 
a k X n matrix A is the unique k x n matrix H of the following form 

/ ... hi /li,2 . . . /ll,n-fc \ 
... /l2 ••• h2,n-k 

\00...0 ... hk J 

where < hij < hi, such that there exists U G G1„(Z) with A = HU, see e.g. |^. We have 
then the following result. 

Theorem 6 The density of k x n matrices whose Hermite normal form is the block matrix 

[Ofcx(n-fc) h] is Ak^n- 

Proof: The result follows from the fact that the set E oi k x n unimodular matrices over Z 
is the equal to the set L of k x n matrices whose HNF is the block matrix [Ofcx(n-fc) ^k]- Let 
us show that E C L. If A & E then there exists a unimodular square matrix B such that 
A consists in the last k rows of B. The HNF of B is the identity matrix In- Indeed, it is a 
upper diagonal square matrix with integer entries whose determinant is 1, which forces the 
diagonal entries to be 1. The size condition of the coefficient of the HNF forces the entries 

6 



above the diagonal to be 0. Thus we have BU^^ = In for some U € G1„(Z), which gives 
AU~^ = [Oj.x(n-k) Ik]- The uniqueness of the HNF shows that A G L. Let us show now 
that L C E hy showing that the gcd of the k x k minors of a matrix in L are coprime. The 
equation A = [Oi^x(n-k) Ik]U shows that the k x k minors of A are equal to the k x k minors 
of [Ofcxfn-ik) ^k] which are or 1, and thus coprime. This finishes the proof of the corollary. 

D 

Taking into account that the Smith normal form of a matrix [3] is obtained from the 
Hermite normal form via row operations, an immediate consequence of this theorem is: 

Corollary 7 The density of k x n matrices whose Smith norm^al form is the block matrix 
[Okyi{n-k) h] is dk,n- 
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